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Abstract 

The Skyrme model is a low energy, effective field theory for QCD which when coupled to 
a gravitational field provides an ideal semi-classical model to describe neutron stars. We use 
the Skyrme crystal solution composed of a lattice of a-like particles as a building block to 
construct minimum energy neutron star configurations, allowing the crystal to be strained 
anisotropically. We find that below 1.49 solar masses the stars' crystal deforms isotropically 
and that above this critical mass, it undergoes anisotropic strain. We then find that the 
maximum mass allowed for a neutron star is 1.90 solar masses, in close agreement with a 
recent observation of the most massive neutron star yet found. The radii of the computed 
solutions also match the experimentally estimated values of approximately 10km. 
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1 Introduction 



Neutron stars are stars that have collapsed under intense self gravitational pressure to the 
point where all electrons are squeezed into nuclei, hence forming a large cluster of neutrons 
with a typical radius of about 10km. A neutron star can thus be seen as a gigantic nuclei 
that is electrically neutral but is strongly affected by the gravitational field that it generates. 
For lack of a unified theory of strong interactions and gravity one has to resort to finding 
an approximate theory that allows us to describe such a system. 

One such theory is the Skyrme model. Originally proposed by Skyrme in 1961 [T], [2] 
as a nonlinear theory of pions to describe strong interactions, it was later shown by Witten 
[3] to be an approximate, low energy, effective field theory for QCD which becomes more 
exact as the number of quark colours becomes large. 

Each solution of the Skyrme model is characterised by an integer valued topological 
charge which can be identified with the baryon number B. The simplest solution, B = 1, 
is made out of a so-called Skyrmion and corresponds to a proton or neutron. At the semi- 
classical level, the Skyrme model does not distinguish between a neutron and a proton. 
Moreover, as the model does not include the electroweak interaction, all Skyrmions are 
electrically neutral. 

The B = 1 solution of the Skyrme model is the only exact stable solution that can be 
computed easily [T]. Solutions with larger values of B can only be computed numerically 
[3], [5] and these solutions have been shown to successfully describe various nuclei and their 
properties [6]. 

Moreover, one can also compute crystal-like solutions made out of an infinite number 
of Skyrmions. In particular, it has been shown that the Skyrme solution with the lowest 
energy per Skyrmion corresponds to a cubic lattice where each lattice unit has a topological 
charge B = 4 [3. These solutions can thus be seen as a crystal of a particles. 

This lattice of Skyrmions looks thus as the best building block to describe a neutron 
star as its has the lowest possible energy per baryon. Yet, one must first estimate if a star 
could instead be made out of a liquid or gas of Skyrmions. The temperature of a neutron 
star, a few years after its creation, cools down to an approximate temperature of about 
lOOeV « 10®/C [g. While this looks fike a very high temperature compared to the binding 
energy of an electron around a nucleus, this energy is quite small from a nuclear point of 
view. Indeed, the lowest excited state of an a particle, for example, is 23.3MeV [9] and the 
lowest vibration mode of a B = 4 Skyrmion is of the order of lOOMeV [10|.[11). Even under 
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intense gravitational energy, Walhout [12] showed tliat the excitation energy of a lattice of 
B — 1 Skyrmions is also of the order of lOOMeV. This points out that the neutron star will 
be in a solid phase rather than a liquid or a gas and that the thermal energy will only excite 
acoustic phonon modes. It is thus natural to model a neutron star as a lattice of B = 4 
Skyrmions. 

Before we proceed we must also question the possibility of having an atmosphere around 
the star, and to estimate its height if it turns out not to be small. At the surface of a 
neutron star twice the mass of the sun, the gravitational acceleration is g ~ 2.6 x lO^^ms^^. 
It is then easy to compute that the average height that an a particle with a thermal energy 
of lOOeV will be able to jump is of the order of 1mm, i.e. much smaller than the radius of 
the star. We can thus consider that such an atmosphere is extremely thin and assume in 
our model that the neutron star is fully made out of a solid. 

Having established that we can use a Skyrme crystal as the building block to describe 
a neutron star we will proceed as follows. First of all, we will use the equations of state 
computed by Castillejo et al. [7] for the _B = 4 crystal when the lattice is deformed asym- 
metrically. Following Walhout [T2] we will then use a Tolman-Oppenheimer-Volkoff (TOV) 
equation [13], [2], generalising it to allow for matter to be anisotropic [15]. The TOV 
equation describes the static equilibrium between matter forces within a solid or fluid and 
the gravitational forces self-generated by the matter for a spherically symmetric body. 

Combining the TOV equation with the equations of state of the Skyrme crystal, we will 
be able to find configurations that are spherically symmetric distributions of anisotropically 
deformed matter in static equilibrium and so are suitable to model neutron stars. Solving 
these equations numerically for large stars we will show that below a critical mass of 1.49 
solar masses (Mq = 1.98892 X 10^"kg) aU neutron/Sky rmion stars are made out of an 
isotropically strained crystal. We will then show that at this critical mass, there is a phase 
transition and that heavier stars are made out of an anisotropically deformed crystal that 
is less strained radially than tangentially. We will also show that these stars can have a 
mass of up to 1.90Mq. Finally, we will investigate the impact of adding a mass term to the 
Skyrme model and describe what happens to a star when its mass is increased above its 
maximum value. 

Using Skyrmions to model neutron stars is not new and has been performed previously 
in several ways. First of all, Walhout used a lattice of _B = 1 Skyrmions [12] to describe 
a neutron star. He then improved his results by considering a lattice of B = 4 Skyrmions 
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|16j . In both cases he assumed an isotropic compression of the lattice, assuming a gas-hke 
phase, and he used numerical solutions of the model to estimate the stress tensor. The 
maximum mass he obtained for the neutron star was 2.57 Mq. Later, Jaikumar and Ouyed 
[17| considered the equation of state for a neutron star based on a Skyrme fluid and obtained 
a maximum mass of 3.6M0 . The main difference between these 2 approaches and ours is that 
they assumed an isotropic fluid of Skyrmions whereas we consider a solid crystal allowed to 
deform anisotropically, i. e. be compressed differently in the radial and tangential directions 
of the star. In our previous papers [TH], [H], we computed minimal energy Skyrmion stars 
made out of layers of 2 dimensional Skyrme lattices. This allowed us to use the rational 
map ansatz [2D] to minimise the energy directly but resulted in relatively small stars with 
a maximum mass of O.574A'f0. This was mainly due to the fact that the field transition 
between the different layers in our ansatz over estimated the energy of the configuration 
and that the energy per baryon in each layer of the ansatz was also larger than that of the 
crystal of B = 4 Skyrmions. 



2 Skyrme Crystals 

The Skyrme model PP, [2] is described by the Lagrangian 

Csk = ^TrCV^C/VC/-^) + J^Tr[(V^[/)C/-\ (V.[/)C/-^]^ + ^4^Tr(f/ - 1), (1) 
to j2e-^ 8 

where here U, the Skyrme field, is an SU{2) matrix and F^, e and are the pion decay 

constant, the Skyrme coupling and the pion mass term respectively. In the Lagrangian ([T]) 

the V are ordinary partial derivatives in the absence of a gravitational field and become 

covarient derivatives when the Skyrme field is coupled to gravity. The results summarised 

in this section all relate to the pure Skyrme model without gravity. 

The Skyrme field is a map from R'^ to S^, the group manifold of SU{2), but finite energy 
considerations imply that the field at spatial infinity should map to the same point, meaning 
the Skyrme field is a map between two three-spheres. Such maps fall into homotopy classes 
indexed by an integer, known as the topological charge, which is interpreted as the baryon 
number, B. The topological soliton solutions, known as Skyrmions, are identified as baryons 
with an a particle described by a _B = 4 Skyrmion solution. 

Here we will be considering the zero pion mass case where — with section 4.3 
describing the effects of its inclusion. The two other Skyrme parameters, F-n and e can be 
obtained in different ways. Skyrme first evaluated them by taking the experimental value 
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of the pion decay constant = 186MeV and then fitting the mass of a Skyrmion to that of 
a proton and obtained e = 4.84. Later Adkins, Nappi and Witten [21] quantised the B = 1 
Skyrmion to fit the parameter values to the mass of the nucleon and the delta excitation and 
obtained F,r = 129MeV and e = 5.45. These later values were the ones used by Castillejo et 
al. [7] to compute the energy of the deformed B — 4 crystal and we will thus use them too. 

The solution of the Skyrme model with the lowest energy per baryon has been shown to 
be a face-centred cubic (fee) lattice of Skyrmions [22], [7]. Each unit cell is a cube of side 
length a with a baryon number of _B = 4 and can therefore be considered as an a particle. In 
the context of a neutron star, we will be able to interpret each B = 4 crystal component as 
being 4 neutrons as the Skyrme model does not distinguish between neutrons and protons. 

Castillejo et al. [7] also investigated the energy of dense Skyrmion crystals where the 
configuration was not a face-centred cubic lattice but rather a lattice where the aspect ratio 
of the unit cell, B = 4 a particle, of side a was altered so that it becomes rectangular with 
aspect ratio r"^. This means that in the x and y directions the lattice size becomes ra and 
in the z direction, a/r^. As in Castillejo et al. we use the measure p = r— 1/r to describe 
the deviation away from the face-centred cubic lattice symmetries which have p = 0. 

The numerical solutions found in [7] provide an equation for the dependence of the 



energy of a single Skyrmion, E{L,p), on its size, L 
number density, and its aspect ratio measure, p. 



-1/3 



where n is the Skyrmion 



E{L,p) = Ep=o{L) + EoHL)p'' + PiL)p^ + -y^p^ + S{L)p' + , 



where the coefficients are given by 

Ep=o{L) = Eo 



5{L) 



0.474 



L 



a{L) = 0.649 - 0.487— + 0.089 ^ 



+ 0.0515 
Lo 



0.300 + 0.006^ - 0.119^, 

Lo L 

-1.64 -f 0.78^ -f 0.71^, 
Lo L 



0.53-0.55 



L 



(2) 

(3) 
(4) 
(5) 
(6) 
(7) 



Here Eo = 727.4MeV and Lo = 1.666 x 10~^^m. The equation can be extended to include 
lower densities [7] but they are not of interest here where we are only considering densities 
higher than the minimal energy crystal. Notice that for any value of L the minimum energy 
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occurs at the face-centred cubic lattice configuration, p = 0, and the global minimum is 
reached for L = Lq. 

3 TOV Equation for Skyrmion Stars 

Using equation ((2)1 relating the energy of a Skyrmion to its size and aspect ratio we will 
now investigate how one can describe a neutron star using a Skyrme crystal and how this 
crystal is deformed under the high gravitational field it experiences. 

In our numerical work we denote Xr as the Skyrmion length in the radial direction of 
the star and At as the Skyrmion length in the tangential direction. These parameters and 
the parameters L and p used in ^ are related as follows 

L = (A,AAt)i andp= (^^y - (^^y . (8) 

To construct a neutron star we consider a spherically symmetric distribution of matter 
in static equilibrium with a stress tensor that is in general locally anisotropic. Spherical 
symmetry demands that the stress tensor, T^', is diagonal and that all the components are 
a function of the radial coordinate only. We denote this stress tensor as 

= diag(p(r),p.(r),pfl(r),p^(r)), (9) 

and consider that, again due to spherical symmetry, pe{r) — Ptt,{r) which we will denote by 
Pt{r) — pe{r) — Pt/>{r). The quantities Pr{r) and pt{r) describe the stresses in the radial and 
tangential directions respectively while the quantity p(r) is the mass density. 

A generalised TOV equation [13], [14] to describe a spherically symmetric star composed 
of anisotropically deformed matter in static equilibrium has been studied previously [13] and 
we summarise it now. 

The metric for the static spherically symmetric distribution of matter can be written in 
Schwarzschild coordinates as 

ds^ = e-^^Ut^ - e^^^^dr^ - r^dO^ - sin^ (10) 

where e"'''' and e^'"^' are functions of the radial coordinate that need to be determined. The 
combination of this metric and the matter distribution, described by the stress tensor ([9]), 
must be a solution of Einstein's equations 

Gab ~ Rab — -RQab = ^TlTab, (H) 



6 



where we have set G = c = 1. After calculating the Ricci tensor and Ricci scalar from the 
metric we find 



-'(^7--)+- = (12) 

(z: + Z2 ] - Z2 = s^P- (13) 

e-^(l."-lAV + i(.')^ + ^^^^) = (14) 



A' 



-A 





1 


= 87rp 




1 




- A' 


;j 


= 87rpt 


2r 





Equation (|12p can be rewritten as 

(re"^)' = 1 - Svr/or^ (15) 

and integrated to give 



(16) 



-A 1 2m 

e =1 

r 

where m = m(r) is defined as the gravitational mass contained within the radius r and can 
be calculated by 

r 2 

m = iirr pdr. (17) 
Jo 

We can now substitute equation p6|l for e^^ into equation (|13p to find 

1 , m + 4-n-r^Pr , , 

r ^ r(r-2m) ' ^^^^ 

The generalised TOV equation that we will use to find suitable neutron star configura- 
tions can now be obtained by differentiating equation (|13|) with respect to r and adding it 
to equation p4p to find 

^=~iP + Pr)'^ + -iPt-Pr) . (19) 

ar Z r 



Now, substituting (|18l) into p9|) . we get 

^ = -ip + Pr)^^^J^^ + %.-pr). (20) 

dr r(r — 2m) r 

For this generalised TOV equation to be solvable two equations of state need to be 
specified, Pr = Pr{p) and pt =pt{p), where, as argued above, we are able to use a zero 
temperature assumption. 

We also need to specify appropriate boundary conditions. First, we must require that 
the solution is regular at the origin and impose that m(r) — !> as r ^> 0. Then Pr must 
be finite at the centre of the star implying that 1/' ^ as r ^ 0. Moreover, the gradient 
dpr/dr must be finite at the origin too and so (pt — Pr) must vanish at least as rapidly as 
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r when r —> 0. This implies that we need to impose the boundary condition pt = Pr at the 
centre of the star. 

The radius of the star, 7?, is determined by the condition Pr{R) = as the radial 
stress for the Skyrmions on the surface of the star will be negligibly small. The equations, 
however, do not impose that pt{R) vanishes at the surface. One should also point out that 
physically relevant solutions will all have Pr,Pt > for r < R. We note that an exterior 
vacuum Schwarzschild metric can always be matched to our metric for the interior of the 
star across the boundary r — R as long as Pr{R) = 0, even though pt{R) and Pr{R) may 
be discontinuous, implying that the star can have a sharp edge, as expected from a solid 
rather than gaseous star. 

As we are considering Skyrmion matter at zero temperature the equations of state that 
will be used in finding suitable neutron star configurations can be calculated from equation 
([2]) which depends on the lattice scale L, and aspect ratio, p, which are both functions of 
the radial distance form the centre of the star, r. From the theory of elasticity we then find 
that the radial and the tangential stresses are related to the energy per Skryrmion, Eq ([2]), 
as follows 

I dE 1 dE 

Using the generalised TOV equation (|20|) and the two equations of state H2ip . a mini- 
mum energy configuration for various values of the total baryon number can be calculated 
numerically. The minimum energy configuration is defined as the minimum value of the 
gravitational mass, Mq, 

l-R 

Ma = m{R) = m(oo) = / '^ivr'^pdr, (22) 



where R is the total radius of the star and 

E 



(23) 



We now need to minimise Ma as a function of A,- and At which both depend on r. 
To achieve this, we will first assume a profile for At(r) and compute Mq for this profile 
as described below. We will then determine the configuration of the neutron star, with a 
specific baryon charge, by minimising Mq over the field At. This can be easily done using 
the simulated annealing algorithm. 

To compute Mq we notice that at the origin, one can use (|2ip to determine Pr(0) and 
pt(0) from the initial values of Ar(0) and At(0). Then the integration steps can be performed 



8 



as follows. Knowing Xr{r) and At(r) one computes p{r) using (|23|l and m(r) using (|17p . 
Then, knowing Pr(''), Pt{r), p(r) and m(r) one can integrate (|20|l by one step to determine 
Pr('' + dj*)- One can then use (|2ip to determine Xrir + dr) and as the profile for At(r) is 
fixed, one can proceed with the next integration step. 

One then integrates (|20p up to the radius R for which Pr{R) ~ 0; this sets the radius of 
the star. In our integration, we used a radial step of 50m. 

One must then evaluate the total baryon charge of the star using 

B = / ^r- , dr (24) 

/ /I 2Gm M /9 ^ ' 

where 



(1-^)^/^ 



- MOW ^^'^ 

and rescale At to restore the baryon number to the desired value. One then repeats the 
integration procedures until the baryons charge reaches the correct value without needing 
any rescaling. 



4 Results 

4.1 Stars Made of Isotropically Deformed Skyrme Crystal 

We found that up to a baryon number of 2.61 x 10^^, equivalent to 1.49Af0, the minimum 
energy configurations are all composed of Skyrme crystals that are isotropically deformed, 
with Ai(r) = Ar(r) across the whole radius of the star. It can be shown that this indeed has 
to be the case as we can prove that if it is possible to find an isotropic Skyrme crystal solution 
then that solution will be the minimum energy configuration. Such isotropic solutions can 
only be found up to a baryon number of 2.61 x 10^^. Corchero [23] used a similar proof for 
a quantum model of neutron stars and we adapt this here for our Skyrme crystal model: 

// there is a locally isotropic, stable solution to the generalised TOV equation H19p with 
mass M and total baryon number N, then all locally anisotropic solutions that have the 
same total baryon number N, in the neighbourhood of that stable solution, will have a mass 
not smaller than M . 

To prove this, we first note that stable solutions for a given baryon number at zero 
temperature, by definition, have a mass that is not greater than that which could be achieved 
by any variation of the density that preserves the baryon number |24j . 
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We then consider two changes to a stable solution with mass Mi, baryon number A'^i, 
density pi(r) and number density ni(r). The first involves changing the density from pi(r) 
to P2{r) while keeping the total baryon number constant and preserving locally isotropy. 
This will result in a configuration that has a mass M2 that is greater than or equal to the 
mass Ml of our initial configuration as that was defined as the minimum mass solution. This 
new configuration will have a different number density, 712 (j*), but the same total baryon 
number, A''i, by assumption. 

The second change involves introducing local anisotropy while the mass, M2, remains 
the same, as does the density, p2{r). In order to keep p2{r) constant when we alter the 
configuration so that it is now made of anisotropic Skyrme crystal the number density 
must also be altered. A change from an isotropic to an anisotropic Skyrme crystal involves 
increasing its energy, and therefore mass, so to keep its mass density constant we need to 
reduce the number of Skyrmions to the new number density n3(r) < 712 (r), meaning the 
total baryon number is now N3 < Ni. 

The two changes described have the effect of firstly increasing M without changing 
and then, secondly, decreasing A'' without altering M. We know that M is a monotonically 
increasing function of A' for isotropic Skyrme crystal stable star configurations, so we have 
proved that moving from isotropic to anisotropic Skyrme crystal configurations increases 
the energy for a given baryon number so does not produce a minimum energy solution. 

3e+57 I 1 1 1 1 1 1 




I 1 1 1 1 1 1 

0.4 0.6 0.8 1 1.2 1.4 1.6 

Skyrmion Length at the Origin, L(r=0) (fm) 

Figure 1: Total baryon number as a function of the size of the Skyrmions at the centre of 
star, L{r = 0). 
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The above proof however does not rule out the existence of anisotropic Skyrme crystal 
solutions for those baryon numbers for which there does not exist an isotropic Skyrme crystal 
solution and such configurations will be discussed in the next section. 

To confirm the results obtained for isotropically deformed crystals, we will now determine 
the properties of these symmetric stars by imposing that symmetry, i.e. pt — p,-. In this 
case the problem simplifies greatly and the TOV equation H2Up reduces to 

= + r{r-2m) ^^O) 
Using this standard TOV equation, a central Skyrmion length \t{r — 0) — Xr{r — 0) = 
L{r — 0) can be specified at the centre of the star. The equation can then be numerically 
integrated over the radius of the star using the Skyrmion energy equation ((2| with 

where, as we are only considering isotropic Skyrme crystal deformations, Xt — Xr and p = 
in the energy equation. This was done using a fourth order Runge Kutta method over 
points every 20m. Notice that this did not require the explicit minimisation of Mq- Figure 
[T] shows a plot of the total baryon number of the star against its Skyrmion length at the 
center, L{r = 0), calculated using this method. 

We found that isotropic Skyrme crystal solutions can be found up to a baryon number 
of 2.61 X 10^*^, which is equivalent to a mass of 1.49Mq. This agrees with the results that 
we found from our minimisation procedure using the generalised TOV equation that allows 
for anisotropic Skyrme crystal deformations. 

Table[T]shows some of the properties of the minimum energy solutions for various baryon 
numbers obtained from the energy minimisation of the generalised TOV equation. The 
results are in perfect agreement with the results obtained by solving the isotropic TOV 
equation (|26|) . The quantity Smin is the minimum value, over the radius of the star, of 

5(r) = e-^^'-^ = l-^, (28) 



which appears in the static, spherically symmetric metric (|10p that we are considering. The 
zeros of S{r) correspond to singularities in the metric, or in other words, to horizons. Had 
Smin been negative, we would have concluded that the neutron star would have collapsed 
into a black hole, but this never occurred. 

We note that the solutions are energetically favourable as the energy per baryon decreases 
when the total baryon number increases, indicating that the solutions are stable. They 
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B 


Total Energy (J) 


Encray/S (J) 


Mass/Mg 


i?(m) 


Smm 


1.0 X 10^^ 


1.16210 X 10^5 


1.16210 X 10-1° 


0.00649160 


2219.20 


0.991503 


1.0 X 10^6 


1.15114 X 10"*^ 


1.15114 X 10^1" 


0.0643083 


4714.35 


0.959976 


2.0 X 10^6 


2.28551 X lO'^s 


1.14276 X 10-^" 


0.127680 


5875.04 


0.936375 


4.0 X 10^6 


4.51669 X 10"*^ 


1.12917 X 10^1" 


0.252325 


7266.13 


0.897929 


6.0 X 10^6 


6.70497 X 10^6 


1.11750 X 10-1° 


0.374573 


8177.42 


0.865580 


8.0 X 10^6 


8.85463 X 10^6 


1.10683 X 10-1° 


0.494664 


8852.67 


0.835232 


1.0 X 10" 


1.09679 X 10-^^ 


1.09679 X 10-1° 


0.612721 


9379.47 


0.808115 


1.2 X 10" 


1.30461 X 10"*^ 


1.08718 X 10-1° 


0.728823 


9798.86 


0.781969 


1.4 X 10" 


1.50899 X 10^^ 


1.07785 X 10-1° 


0.842997 


10133.2 


0.755523 


1.6 X 10" 


1.70994 X 10^^ 


1.06871 X 10-1° 


0.955258 


10394.6 


0.730148 


1.8 X 10" 


1.90741 X lO'^^ 


1.05967 X 10-1° 


1.065578 


10588.7 


0.704181 


2.0 X 10" 


2.10132 X lO'^^ 


1.05066 X 10-1° 


1.173903 


10714.6 


0.677181 


2.2 X 10" 


2.29147 X 10'^^ 


1.04158 X 10-1° 


1.280129 


10761.8 


0.649383 


2.4 X 10" 


2.47750 X 10^^ 


1.032293 X 10-1° 


1.38406 


10694.5 


0.619124 


2.6 X 10" 


2.65860 X 10^^ 


1.022536 X 10-1" 


1.48522 


10367.5 


0.577658 



Table 1: Properties of the isotropic minimum energy neutron star configurations for various 
baryon numbers. 
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correspond to the solutions to the right of the maximum in figure [T] with sohitions to the 
left being unstable with a higher energy per baryon for a given baryon number, and therefore 
not found by the energy minimisation procedure. 

The neutron star solutions which have masses larger than the mass of the Sun have 
radii of about 10km, which very much matches the experimental estimates of the radii of 
observed neutrons stars. Notice also that the largest neutron star, in our model, has a mass 
of approximately 1.28Af0, and above that value, the radius of the stars decreases with their 
mass (see table [T] and figure (2]). 




Mass (Mq) 

Figure 2: Radius of the neutron star solutions as a function of their mass (solid line), and that 
of the maximum mass solution (cross). 

We now consider the structures of these isotropic Skyrme crystal stars, in particular we 
consider the case of a star with a mass of I.4OM0, a typical mass for a realistic neutron 
star, equivalent to a baryon number of 2.44 x 10^^, although all the isotropic Skyrme crystal 
minimum energy solutions show the same qualitative behaviour. Figure [3] shows the size of 
the Skyrmions, L{r), over the radius of the star. As expected the Skyrmions are deformed 
more towards the centre of the star than at the edge, increasing the Skyrmion mass density 
by a factor of 4.44. Due to this decrease in the size of the Skyrmions as we reach the centre 
of the star the stress is higher at the centre and decreases towards zero at the edge of the 
star as imposed by the boundary conditions. 

The isotropic Skyrme crystal solutions have a Smin that is always greater than zero so 
the configurations do not collapse into black holes. Figure O also shows how the value of 
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Radius, r (km) 

Skyrmion Length, L(r) S(r) 

Figure 3: Variation of the size of the isotropic Skyrmions, L{r), (solid line) and of the metric 
function S{r) (dotted line) over the radius of a star of mass 1.4OM0. 



S{r) varies over the radius of the star. 

4.2 Stars Made of Anisotropically Deformed Skyrme Crystal 

Having shown in the previous section that no isotropic Skyrme crystal solutions exist for 
baryon numbers larger than 2.61 x 10^^, we will now show that anisotropic solutions do 
exist. 

Table [2] shows some of the properties of the anisotropic minimum energy Skyrme crystal 
solutions for various baryon numbers obtained using the generalised TOV equation. We 
found solutions in this way up to a baryon number of 3.25 x 10*^, corresponding to I.SIMq, 
after which the numerical energy minimisation procedure became difficult to implement. 
However by using a similar simulated annealing process to maximise the baryon number, 
rather than minimise the energy for a particular baryon number, we found anisotropic 
Skyrme crystal solutions up to a baryon number of 3.41 x 10^^, equivalent to 1.9OAf0. 
At this maximum baryon number solution there is only one possible configuration of the 
Skyrmions, as any modification to it results in a decrease in the baryon number, hence it is 
the minimum energy solution. Above this baryon number, solutions do not exist. 

As in the case of isotropic Skyrme crystal deformations we find that the solutions are 
energetically favourable as the energy per baryon decreases as the total baryon number 
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B 


Total Energy (J) 


Energy/B (J) 


Mass/Afg 


i?(m) 




2.65 X 10" 


2.70277 X 10^7 


1.01991 X 10-1° 


1.50990 


10091.8 


0.559060 


2.70 X 10" 


2.74605 X 10''^ 


1.01706 X 10-1° 


1.53408 


9555.51 


0.526465 


2.75 X 10" 


2.78943 X 10'^'' 


1.01434 X 10-1° 


1.55832 


9460.46 


0.514207 


2.80 X 10" 


2.83310 X 10*7 


1.01182 X 10-1° 


1.58271 


9456.89 


0.506402 


2.85 X 10" 


2.87706 X 10^7 


1.00949 X 10-1° 


1.60727 


9456.46 


0.498735 


2.90 X 10" 


2.92133 X 10^7 


1.00735 X 10-1° 


1.63200 


9457.92 


0.491152 


2.95 X 10" 


2.96592 X 10^^ 


1.00540 X 10-1° 


1.65691 


9460.65 


0.483633 


3.00 X 10" 


3.01087 X 10^7 


1.00362 X 10-1° 


1.68202 


9465.06 


0.476231 


3.05 X 10" 


3.05619 X 10*7 


1.00203 X 10-1° 


1.70734 


9469.97 


0.468880 


3.10 X 10" 


3.10191 X 10^^ 


1.00062 X 10-1° 


1.73288 


9475.76 


0.461631 


3.15 X 10" 


3.14807 X 10^7 


9.99388 X 10-11 


1.75867 


9481.95 


0.454438 


3.20 X 10" 


3.19472 X 1047 


9.98351 X 10-11 


1.78473 


9489.04 


0.447382 


3.25 X 10" 


3.24191 X 10^7 


9.97510 X 10-11 


1.81109 


9496.62 


0.440435 



Table 2; Properties of the anisotropic minimum energy neutron star configurations for various 
baryon numbers. 
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increases, indicating stable solutions. As the baryon number is increased towards its max- 
imum value of 3.41 x 10*^ the energy per baryon begins to level off and we find that the 
ma:ximum baryon number has the lowest energy per baryon, as in the isotropic case. 

We can see that the configurations we have constructed do not collapse into a black hole 
by noticing that the values of Smin are always positive, as shown in figure [l] 




Mass (Mq) 

Figure 4: Smin of the neutron star solutions as a function of their mass. The maximum mass 
solution is shown as a cross. 

Figure [2] shows a plot of the mass radius curve for both the isotropic and anisotropic 
Skyrme crystal cases, with the mass in units of AIq. As stated above, large isotropic crystal 
neutron stars have a radius that decreases as the mass increases. We can clearly see in figure 
[21 that at the critical mass of 1.49M0, the radius keeps decreasing as the mass of the star 
increases. Moreover, we also observe a sharp drop of radius just over I.SMq followed by a 
plateau at about 9.5km. 

By considering anisotropic as well as isotropic Skyrme crystal solutions we have extended 
the mass range over which solutions can be found, finding masses up to 28% above the 
maximum mass of the isotropic case. This is an interesting finding because isotropy of matter 
is often taken as an assumption when studying neutron star models, including the Skyrme 
crystal case considered in [12|.[16|.[17j and a maximum mass is then derived. We have shown 
that by not assuming isotropy and instead allowing anisotropic matter configurations the 
maximum mass can be increased by a significant amount. In this simple Skyrme crystal 
model the maximum mass found is equivalent to l.QOAf© and the recent discovery of a 
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(c) (d) 

Figure 5: Skyrmion lengths A, (r) (solid line), At(r) (dashed line) and L{r) (dotted line) for a) 
Largest neutron star {R ~ lO.Skm): AI — 1.28A'/q b) Heaviest isotropic neutron star: A/ = 
1.49Mq (all lengths coincide as they are made of isotropically deformed crystal); c) Densest 
neutron star: M = 1.54Mq; d) Heaviest neutron star: M = 1.90Mq. 
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1.97±0.04Mq neutron star [25], the highest neutron star mass ever determined, makes this 
an encouraging finding, especially when we consider that including the effects of rotation 
into our model will increase the maximum mass found, by up to 2% for a star with a typical 
3.15ms spin period pS] . 

Figure [5] shows a selection of plots of the Skyrmion lengths \r and At and the Skyrmion 
size L, equation ([SI, over the radius of the star for four special stars: the largest star, 
with radius R = 10.8km and mass M = 1.28M0 (figure I5a|) : the heaviest isotropically 
deformed star AI = 1A9Mq (figure [5b)) : the densest neutron star, M — 1.54M0 (figure 
[5c)l and the heaviest neutron star, M — 1.90Mq (figure [5d[) . The first two are made out 
of an isotropically deformed crystal, while the last two are anisotropically deformed and 
one notices that the amount of anisotropy increases as the mass increases (the divergence 
between \r and At increases). Throughout this paper, we will use these four special stars 
as examples to illustrate various properties of the neutron stars. 

As the maximum mass is approached the gradient of the profile of tangential Skyrmion 
lengths over the radius of the star becomes smaller and we note that physically mean- 
ingful stars composed of anisotropically deformed crystal should have dXt/dr > |27[ . 
This confirms that the minimum energy solution for the maximum mass found, I.9OM0, 
for anisotropic Skyrme crystal solutions is the configuration with a constant tangential 
Skyrmion length as illustrated in figure [5dl 

The generalised TOV equation imposes that the sizes of the Skyrmions are equal in 
all directions at the centre of the star, but away from the centre, for all the anisotropic 
Skyrme crystal solutions, we find that the amount of Skyrmion anisotropy increases as we 
move towards the edge of the star, reaching the maximum at the edge. The Skyrmions are 
deformed to a greater extent in the tangential direction in agreement with the value of the 
aspect ratio, p, being negative over the values where Ar 7^ At- 

As expected, the profiles for Ar and At show that the mass density at the centre of the 
star is higher than at the edge, decreasing monotonically as the radial distance increases. 
This is shown by figure [B] for the largest, heaviest isotropic, densest and maximum mass 
solutions. 

In figure [7| one can see how the lengths of the Skyrme crystal Ar and At vary with 
the mass of the star both at the center (r = 0) and the edge of the star (r = R). For 
isotropically deformed stars, Xr{R) = At(-R) is constant and corresponds to the minimum 
energy Skyrme crystal in the absence of gravity. Not surprisingly, Xr{0) = At(0) decreases 
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Figure 6: Mass density p{r) for: a) Largest neutron star (i? = 10.8km): M ~ 1.28Mo (solid 
line) b) Heaviest isotropic neutron star: M = 1.49M0 (dashed line); c) Densest neutron star: 
M = 1.54Mo (dotted line); d) Heaviest neutron star: M — LOOM© (dash dotted line). 

steadily as the mass of the star increases, showing that the density at the center of the star 
increases. Once the phase transition has taken place and the star is too heavy to remain 
isotropically deformed, we observe that Ar(0) — At(0) drops sharply to a local minimum, 
reached for M ~ 1.54Mq. Meanwhile, \r{R) and Xt{R) remain nearly identical. Beyond 
the minimum of Ar,t(0), Xr{R) and Xt{R) start to diverge sharply; Xr{R) decreases slightly 
in value while Xt{R) decreases rapidly. These stars are thus much more compressed in the 
tangential direction than in the radial one. As seen on figure [5dl Xt{R) = At(0) for the 
maximum mass neutron star. 

Another property of a neutron star worth considering is the speed of sound. To compute 
it one needs to know how the energy of the crystal varies when it is deformed in the direction 
of wave propagation. Using ([2]) we can thus compute the speed of sound in the z direction. 
To compute the speed of sound in the x and y directions when the crystal is deformed we 
need to know how the energy of the crystal varies when the crystal is deformed in all three 
directions independently, an expression we do not have. 

We are thus only able to compute the radial speed of sound inside a neutron star and it 
is given by 




(29) 
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0.8 

1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 

Mass (Mo) 
?4(R) A.r(0)=?4(0) 

Figure 7: Skyrmion lengths at the edge of the star, Xr{R) (solid line) and At(i?) (dashed line), 
and at the center of the star, Xr(0) = At(0) (dotted line), as a function of the star mass. 



where both pr and p are functions of A,- and At given respectively by (|21[) and H23p . Obvi- 
ously, when the crystal inside the star is isotropically deformed, the speed of sound is the 
same in all 3 directions. 

First of all it is interesting to notice that the speed of sound in the minimum energy 
Skyrme crystal, in the absence of a gravitational field, is amazingly large: v — 0.57 c. This 
is the speed of sound at the surface of a neutron star when it is deformed isotropically. 
From figure [H] one sees that Vr increases as one moves towards the center of the star. As Vr 
is directly related to the density of the star, it is not surprising to find that the maximum 
radial speed, Vr — 0.78c, is reached at the center of the densest neutron star, i.e. the one 
with M — 1.54iV/0. As expected, Vr < c everywhere. 

Figure shows how the value of S{r) varies over the radius of the star for, again, the 
largest, heaviest isotropic, densest and maximum mass solutions, showing how the metric is 
altered as r varies. The minimum value of S{r) is always located at the edge of the star, i.e. 
Smin = S{R), and it is presented in figure |4] as a function of the star masses. One sees that 
Smin decreases monotonically as the mass increases, and exhibits a sharp decrease just over 
1.5M0, I.e. just above the critical mass. However S,nin always remains positive, indicating 
that no black hole is formed. 

Figure [10] shows how the total baryon number and the mass of all the solutions found 
are related. As the baryon number increases the effects of gravitational attraction increase. 
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Figure 8: Radial speed of sound, Vr{r) for a) Largest neutron star (R = 10.8km): M = 1.28Mo 
(solid line) b) Heaviest isotropic neutron star: M = 1.49Afo (dashed line); c) Densest neutron 
star: M = 1.54M0 (dotted line); d) Heaviest neutron star: M = 1.90 A/© (dash dotted line). 



0.9 



0.8 



^ 0.7 



0.6 



0.5 



0.4 



4 6 
Radius (km) 



l\/l=1 .28Ma, 
IV1 = 1 .49Mn 
IV1=1 .54Ma 
IV1=1 .90l\/15 



10 



Figure 9: The function S{r) for: a) Largest neutron star {R ~ 10.8km): M = 1.28AfQ (solid 
line) b) Heaviest isotropic neutron star: M = 1.49 A/© (dashed line); e) Densest neutron star: 
M ~ 1.54M0 (dotted line); d) Heaviest neutron star: AI = 1.9OAf0 (dash dotted line). 
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Figure 10: Mass of the neutron star solutions as a function of their baryon number. The maxi- 
mum mass solution is shown as a cross. 

resulting in a slightly lower gravitational mass per baryon than expected from a linear 
relation. 

We note that the minimum value of the aspect ratio, p, for the minimum energy con- 
figurations found is —0.283 and the minimum value of L is 8.11 x 10^^^, both of which are 
within the valid range of values for equation ([5} [7]. 

4.3 Inclusion of the Pion Mass 

Throughout the work described we have assumed a zero pion mass. The inclusion of a 
non-zero pion mass can be considered by including the pion mass term, 




(30) 



in the static Skyrme Lagrangian ((T|, where U is the Skyrme field, F-n is the pion decay 
constant and is the pion mass. Using the cubic lattice of a-like Skyrmions that has been 
considered above one finds that Tt{U — 1) = —2, meaning that the energy arising from 



the pion mass term reduces to 




(31) 



an energy term proportional to the volume of the Skyrmions. 

It can be seen in figure [TT] that including a pion mass of m = 138MeV decreases the 
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Figure 11: Mass of the star as a function of the size of the Skyrmions at the centre, Lq, for zero 
pion mass (solid line) and m ~ 138MeV (dashed line). 



maximum mass of the star by a very small amount from 1.49 to 1A7Mq while also slightly 
decreasing the central density at which this occurs. 

Including a pion mass of m = 138MeV in the simulated annealing process used to 
find the maximum baryon number for the anisotropic Skyrme crystal solutions results in a 
maximum baryon number of 3.34 x 10^^, equivalent to I.88M0, a decrease of O.O2M0 from 
the maximum mass found in the case without a pion mass. 

This gives an indication as to how the pion mass affects the structures of the neutron 
star configurations that can be constructed, and a similar reduction in the maximum mass 
is expected for all the anisotropic crystal solutions, however when the pion mass is included 
it also has the effect of driving the Skyrme crystal lattice away from the half-Skyrmion sym- 
metry [7]. This will be a small effect for the dense Skyrme crystals that we are considering 
because while the pion mass term is the dominant term in the Lagrangian far away from the 
centres of the Skyrmions when they are well separated, in the dense Skyrme crystal there is 
no space away from the centres of the Skyrmions so it becomes less important in affecting 
the field distributions. Its effect will be to reduce the pion mass term, Eq. H3ip . by a small 
amount. 
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4.4 Stars above the Mciximum Mass 

As in other studies of neutron starts based on the Skyrme model, we found a critical mass 
above which solutions do not exist. In other words, when the star is too massive, the 
crystal of which it is made is not capable of counterbalancing the gravitation pull and the 
star then collapses into a black hole. This is indeed what we observed when trying to 
construct solutions above the critical mass: the energy of the configuration kept decreasing 
as the radius of the star decreased and the Smin function became negative, indicating the 
formation of an horizon, and hence a black hole. 

Throughout this work we have assumed a spherically symmetric metric and stress tensor, 
however, these assumptions could be removed and it may be that higher mass solutions could 
be found. We could instead consider an axially symmetric metric, the most general form 
pg] being 

ds^ = a^(V + dz'^) + P^dfjP - j^dt'^, (32) 
when written in cylindrical coordinates. The stress tensor, 

T!^ = diag(p,pi,p2,P3), (33) 

could then be completely anisotropic with pi P2 Ps- Minimum energy solutions to 
Einstein's equations for such a metric and stress tensor could be found by direct minimisation 
of the action of the Skyrme model coupled to gravity or by using an, as yet undetermined, 
axisymmetric form of the TOV equation. Another approach to investigate such solutions 
would be to perturb the spherically symmetric solutions that we have found. Following 
the procedure for doing so described in [28] the exterior metric for an axially symmetric 
solution can be written in Schwarzschild coordinates and, after comparing the exterior 
spherically symmetric Schwarzschild solution to our solutions for the interior metric of the 
star and finding the substitutions necessary to move from one to the other, we can make 
to same substitutions to the axially symmetric exterior metric. This allows us to then 
describe approximately both the metric and the stress energy tensor of the axially symmetric 
solution. To carry out such investigations into axially symmetric static configurations an 
equation analogous to ^ which would relate the energy of the Skyrme crystal to its size 
and deformation in all three directions independently would need to be considered. 

We have also assumed that the stress tensor, = dia,g{p, pr , pe , p^,) , is diagonal, how- 
ever, if shear strains are included in our model off diagonal components would have to be 
introduced. This would also remove the assumption of spherical symmetry altering the 
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configurations found. 

Spfiericai symmetry aiso needs to be removed to consider rotating stars. Tliis wiii resuit 
in configurations above the maximum mass found in this work, by up to 2% for a star with 
a typical 3.15ms spin period [26) . and as neutron stars are known to be rotating, this is an 
important effect to consider. 

5 Conclusions 

Neutron stars are large bodies of matter where the electrons, instead of circling atoms, are 
forced to merge with the nuclei, resulting in extremely dense stars made entirely of neutrons. 
Their temperature, from a nuclear point of view, is very low and this means nuclear matter 
must be considered as a solid rather than a fluid. Moreover, the gravitational pull of the 
star is so strong that the "atmospheric" fluid one might expect at the surface is of negligible 
height. 

In this context, the Skyrme model, known to be a low energy effective field theory for 
QCD [3], is an ideal candidate to describe neutron stars once the model is coupled to gravity. 
The minimum energy configuration of large numbers of Skyrmions is a cubic crystal made of 
B = 4 Skyrmions which correspond to a crystal of a-like particles. We have thus used these 
solutions as a building block to describe the neutron star by combining the deformation 
energy computed in [7] and a generalised version of the TOV equation [13] , [14] , [15] which 
describes the static equilibrium between matter forces, within a solid or fluid, and the 
gravitational forces self-generated by the matter for a spherically symmetric body. 

The key feature of our approach to the problem was to consider the star as a solid 
that could potentially deform itself anisotropically. We then found that below 1.49Mq, 
all stars were made of a crystal deformed isotropically, i.e. the radial strain was identical 
to the tangential one. Above that critical value, the neutron star undergoes a critical 
phase transition and the lattice of Skyrmions compresses anisotropically: the Skyrmions 
are more compressed tangentially than radially. Stars were shown to exist up to a critical 
mass of 1.90Mq, a result that closely matches the recent discovery of Demorest et al. [53] 
who measured the mass of the heaviest neutron star found to date, PSR J1614-2230, to be 
1.97Mq . We also observed that the maximum radius for a Skyrmion star was approximately 
11km, a flgure that matches well the experimental estimations. 

In our model we did not consider the rotational energy of the star which is approximated 
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at about 2% of its total energy. If we included that extra energy, our upper bound would 
thus just fit above the mass of PSR J1614-2230. 

Finally we have also shown that if the mass of a neutron star was to be raised to cross 
the critical mass threshold, it would collapse into a black hole. 
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